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Abstract 

We investigate finite temperature effects in O'Raifeartaigh models with global sym- 
metries which exhibit supersymmetry breaking at a meta-stable vacuum accompanied by 
U{l)ii breaking. The pseudo moduli field is stabihzed at one-loop order at zero temper- 
ature within some coupling parameter region. We analyze the behavior of the parameter 
space according to non-zero temperatures and find that the parameter region which allows 
J7(l)/j breaking is considerably extended at sufficiently low temperature, even though it 
shrinks down at high temperature as expected. We also discuss the thermal history of the 
meta-stable supersymmetry breaking vacuum. 
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1 Introduction 



No one doubts that supersymmetry breaking has been the long-standing problem in high energy 
particle physics and the problem is closely related to the existence of global R-symmetry. The 
relation between supersymmetry breaking and global R-symmetry has been studied after its 
connection was established in [1]. The authors of [1] showed that in order to break supersym- 
metry at a global minimum of the potential in generic models, there must exist global U{1)r 
symmetry. However the f/(l)_R symmetry should be broken at the supersymmetry breaking 
vacua otherwise gauginos can not have non-zero masses that are phenomenologically required. 

A possible solution to this puzzle is to consider a local but not global minimum of the 
potential to break supersymmetry. In |2], Intriligator, Seiberg and Shih studied low-energy ef- 
fective potential of an A/" = 1 supersymmetric QCD (ISS model) which exhibits supersymmetry 
breaking at a local vacuum. This vacuum is meta-stable but can be long-lived, compared to 
the life time of the universe. 

Other models with supersymmetry breaking at local vacua have been studied. It has been 
shown that O'Raifeartaigh models [3] with fields that have R-charges different from or 2 can 
have a meta-stable supersymmetry breaking vacuum without U{l)fj symmetry In these 
models, supersymmetry is broken at tree-level while the symmetry is broken at one-loop 

order. Models with tree-level R-symmetry breaking were investigated in |5]. 

It is useful to generalize these models to ones with global symmetries [6]. Such an exten- 
sion is important since if global symmetries are gauged as gauge groups of the Standard Model, 
supersymmetry breaking is mediated to the visible sector through direct gauge mediation mech- 
anisms 0. These observations open up more possibilities for constructing phenomenologically 
viable models (for example see [8]) and the realization of long-lived local vacua in field theories 
provides a good example to study the string landscape problem [9]. 

It is also interesting to study the thermal effects of these local vacua. The thermal effects 
generically modify the structure of the vacua and change the view of the landscape. Since we 
believe that supersymmetry is broken at early universe after the inflation, it is useful to analyze 
the thermal history of the hidden sector. The thermal history of supersymmetric models with 
supersymmetry breaking has been investigated in the context of ISS model [101 [HI [121 [13] and 
the O'Raifeartaigh type model without U{1)r symmetry [T3]. The thermal effects of a broad 
class of models of gauge mediation were also considered [15]. 

In this paper we study effects of finite temperature on the O'Raifeartaigh type model with 
global symmetry considered in [B] . Classically this model has a fiat direction of local extrema 
and a supersymmetric runaway vacuum. One loop corrections stabilize the flat direction, lead- 
ing to a spontaneous t/(l)_R breaking vacuum in some regions of parameters in the theory. We 

^ However, ISS model and models as in [6] lead to suppressed gaugino mass which is phenomenologically 
unacceptable [5]. Deformed ISS models have been considered in, for instance, [7|. 
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take account of finite temperature effects and show that such parameter regions shrink down 
as temperature increases. We also show that parameter regions which are not allowed at zero 
temperature cause U{1)r breaking at low temperature. We also address the finite temperature 
effects along the runaway direction and discuss the thermal evolution of a local vacuum. 

The organization of this paper is as follows. In the next section, we introduce the model and 
discuss the supersymmetry breaking vacuum and the runaway direction at zero temperature. In 
section 3, we study finite temperature effects of the model. Especially we investigate behavior 
of the parameter region according to temperatures. In section 4, we discuss the thermal history 
of the local vacuum focusing on a specific path to the runaway direction. Section 5 is conclusion 
and discussions. In appendix A, analytic expressions of the effective potential in high and low 
temperature expansions are shown. In appendix B, we give explicit expression of general mass 
matrices. 



2 The model and supersymmetry breaking 

In this section, we provide a brief overview of the generalized O'Raifeartaigh model that exhibits 
supersymmetry breaking without U (1) r symmetry [1]. We then introduce the model with global 
symmetry and show the result obtained in [6] at zero temperature. 
Consider O'Raifeartaigh models with the following superpotential 

W = fX + ^{M'^ + XN'^)¥<I>^, (1) 
det M 7^ 0, det(M + XN) ^ 0, (2) 

where M, N are complex symmetric matrices, / is a real parameter and X, are A/" = 1 chiral 
superfields. We assume that the Kahler potential for $*, X is canonical. The condition of U{1)r 
charge for each field is R{X) = 2 and the other fields satisfy 

R{^') + /?($■'■) = 2 for M'^ 0, 
R{¥) + R{<!>^) = for N'^ ^ 0. 

The supersymmetric vacuum conditions are given by 

= / + ^N'^^'cfP, = {M'^ + N'^X)4P. (4) 

Here 0* and X are the lowest component of the superfields and X. Due to the condition ([2]), 
the equations in dlj can not be solved simultaneously and the model breaks supersymmetry at 
an extremum of the potential: 

0* = 0, X : arbitrary. (5) 



2 



This solution leaves a flat direction for the X field, which is called pseudo modulus. In order 
that this extremum is a local vacuum, there must not exist tachyonic directions. The bosonic 
and fermionic mass matrices at this extremum are given by 

^2 ^ ( {M^ + XN^){M + XN) fN^ \ 

\ fN {M + XN){M^X + N^) y' 

^2 ^ ( {M^ + XN^){M + XN) \ 

"^■^ y [M + XN){M^X + N^) J' ^' 

If / = 0, the point ^ becomes a supersymmetric global vacuum. In this case, all the 0- 
directions are not tachyonic and the eigenvalues of m\ = m?p are positive. On the other hand, 
once we turn on / 7^ terms in the off-diagonal parts in m^, these terms can be considered 
as perturbations to the positive definite diagonal parts. If the off-diagonal parts are smaller 
than the diagonal parts, the eigenvalues of m\ can stay in their positive values. Especially if 
we focus on the small X region, the condition that there are no tachyonic directions is given 
by \fN\ <^ |Mp which was employed in [6]. In addition to a local supersymmetry breaking 
minimum, in general, there is also a supersymmetric vacuum in a runaway direction. The local 
minimum is thus a meta-stable vacuum. 

Since the fields 0* have mass m*-' ~ M^^ + XN^^ at the origin of 0', given the energy 
scale lower than this mass, 0* are integrated out and the effective potential for X is obtained. 
The one-loop effective potential of the pseudo modulus X is given by the Coleman- Weinberg 
potential 



2 2 ■ 

4 1 '^B 4 1 
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where mBjiriF are bosonic and fermionic mass matrices and A is the dynamical cutoff scale. 
The total effective potential Ves{X) is just the sum of Vtrec and where Vtree is the tree-level 
potential. The one-loop effective potential generically fixes the modulus X. 

It was shown that the model ([T]) can dynamically break the f/(l)R symmetry at the super- 
symmetry breaking vacuum by the one-loop corrections if there is a field that has R-charge 
different from or 2 [1]. Once the parameters M, and / are chosen appropriately, the pseudo 
modulus X is stabilized a.t X ^ giving a U {1)r breaking vacuum as X has non-zero R-charge. 

The model ([T]) has been generalized to ones with global symmetry and with more pseudo 
moduli [B]. In this paper, we focus on the simplest model that has global symmetry U{N) and 
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one modulus X which is discussed in [B]. The interaction matrices are given by 
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(9) 



The superpotential is given by 
W - 



+Mj(f)1,)^^_5)a + M50^5)0(_3)a + M30^3)0(_l)a, 



(10) 



where (p'^j-^ , ^ 



1, 



A^) are fundamental, 



-5)a, (P{-3)a, ' 



-l)a 



are anti-fundamental 



representations of U{N) global symmetry and X is a U{N) singlet. The field has an R- 
charge r. We have chosen the real parameters /, N^jN^, M3, M5 and M7 to be positive. 

As we have discussed, this model breaks supersymmetry at the vacuum = for all 
r leaving the flat direction parametrized by X. In addition to this local minimum, there is 
a runaway direction to the true supersymmetric vacuum, which is obtained from the vacuum 
conditions dl]). It is given by 



N3 cps 

N.M, (hi 



(3) 
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(11) 



where the constants (/, A''3, N^, (j)^, ())_^) satisfy the following conditions. 



f + 2N, 



0, 



3(P3(, 



(12) 



The runaway direction is characterized by a parameter e — )■ 0. 

For later convenience, we re-derive the result obtained in [B] for the allowed region of ?7(1)_r 
breaking vacua in the parameter space of /, M and at zero temperature. Because the 
effective potential of the model with the U{N) global symmetry is just the copy of the one 
for the model with iV = 1, it is enough to consider = 1 case to see the one-loop effective 
potential for X. We investigate the allowed regions for [M^/M^, Mj/M^,) with A'3 = A^s = 1 
and f /M^ = 0.001. The last value is chosen to satisfy the condition \M~'^fN\ ^ 1 to avoid 
tachyonic directions for small X. The result is given in Fig. [T]-(a). 

In the next section, we consider the one-loop effective potential at finite temperature and 
study the behavior of the parameter region that allows U{1)r breaking vacua. We will see how 
this parameter region changes according to non-zero temperature. 
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3 Finite temperature effects 



The finite temperature contribution in one-loop effective potential is given by [T7] 



'4 



(13) 



'4 



rj'Vc) = -2X—MmU4>c)/T% 



(14) 



where the first is real scalar and the second is fermionic contributions. Here T is temperature, 
0c are VEVs at vacua and the fermionic degrees of freedom A are given by A = 1 (Weyl) or A = 2 
(Dirac). The functions Jb and Jp are the so-called bosonic and fermionic thermal functions 
given by 



In the following, we consider the effective potential (fT7|) of the model (fTOj) . As discussed in 
the previous section, since the model has the U{N) global symmetry, the effective potential 
f lT7|) is just copy of the one with = 1 case. We therefore consider f|T7|) with A^ = 1 case. 
Since the high and low temperature expansions are useful to see the behaviour of the potential 
analytically and to check the full numerical result, these expansions are collected in appendix 
A. 

For all the results in the following, the effective potential f ll7p is directly calculated by 
using Coleman- Weinberg potential ([8]) and the thermal potentials ( 1T3|) and ( IT^ without any 
approximation formulas. 

We show the parameter region which allows non-supersymmetric U{1)r breaking local min- 
imum at temperatures ranging from T = to T = 1.0 in Fig. [T] In this analysis, we study the 
allowed region against the parameters (M3/M5, M7/M5) near the origin of X, \X\ < 3. Their 
allowed region becomes smaller and the parameters [M^/M^, Mj/M^) have larger values, as 
temperature increases. 

The numerical results in Fig. [T] enable us to see the structure of phase transitions of a given 
vacuum during the cooling process. It is seen that at high temperature the U{1)r symmetry 
is generically recovered, which is expected from the analytic calculations in Appendix A. As 
temperature decreases, with an appropriate choice of parameters, a vacuum that preserves the 
U{l)ji symmetry at high temperature becomes a f/(l)_R breaking vacuum. 




(15) 



(16) 



where /3 = T ^. The one-loop effective potential is therefore given by 



(17) 
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Figure 1: Contour plots of the allowed region of the f/(l)R breaking against the parameters 
{M'i/M^^MT/M^) for |X| < 3 with temperatures ranging from T = to 1.0. We choose values 
of the other parameters as N^, = = 1, //M| = 0.001. 

We note that the numerical result in Fig. [T] has a different behavior from the result of the 
model in [H] where the finite temperature effects of the same type of model as ([T]) without 
global symmetry have been studied. In [H], the contour plots similar to Fig. [T] is shown 
and it shows that once certain parameters come into an allowed region, they remain in the 
allowed region until zero temperature. However, in our model, some parameters entered into 
the allowed region , then it goes out the region as temperature decreases. This fact means that 
there are several phase transitions during the cooling process. For instance, for the case with 
iVg = iVg = 1 M3/M5 = M7/M5 = 0.4 and //M| = 0.001, the U{1)r symmetry vacuum takes 
place near the origin around T ~ 0.25 and it disappears around T ~ 0.16. 

Eventually, this parameter choice realizes the U{1)r symmetry vacuum at zero temperature, 
but at lower temperatures, the evolution of vacuum structure is a bit nontrivial. Next, we will 
show detailed evolutions of the potential at lower temperature. 

Fig. shows the plots of the potentials for X with = = 1, M3/M5 = M7/M5 = 0.4 
and f/M^ = 0.001. As we have discussed, there exists a U{1)r symmetry breaking vacuum 
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(d) T = 0.02 (e) T = 0.015 (f) T = 0.01 



Figure 2: Plots of potentials for X with temperatures ranging from T = 0.2 to T = 0.01. We 
choose Ns = N5 = 1, = M-j/M^ = 0.4 and //M| = 0.001. In the plots the numerical 

value is evaluated with M5 = 1. 

near the origin for this parameter choice at the temperature around T ~ 0.25 and the vacuum 
disappears at a critical temperature as temperature decreases down to T ~ 0.16. This behavior 
can be seen in Fig. [2]-(a), (b) and (c). Although this U{l)Ft. breaking vacuum around the origin 
disappeared at the critical temperature, the U{l)ji breaking minimum re-appears at lower 
temperatures. For instance, at T = 0.16 we find that there exists a U{1)r symmetry breaking 
vacuum around X ~ 13 which is far from the plot region of Fig. [2l As temperature decreases 
from T = 0.16, this vacuum approaches to the origin as can be seen around X ~ 4 in Fig. 
[2l-(d). As temperature further decreases, the t/(l)_R symmetry breaking vacuum develops near 
the origin again (see Fig[2]-(e)) because the Coleman- Weinberg potential becomes dominant 
and the extra vacuum far away from the origin disappears as can be seen in Fig. 4-(f). 

4 Runaway direction in finite temperature 

In the previous section we studied the allowed region of the parameter space where U{1)r 
symmetry is broken. One would be afraid about that thermal effects on the effective potential 
may prevent the state from staying in the vicinity of the meta-stable vacuum during the cooling 
process and the system may escape to the runaway vacuum through the 0-directions. In this 
section we study finite temperature effects on a potential barrier that intermediates between 
the meta-stable vacuum and a runaway direction. 
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It is not easy to make a thorough investigation how the structure of the effective potential, 
which depends not only on X but also is affected by the change of temperature. Instead, 
according to [lOl [H] , we choose a particular path which passes from the meta-stable vacuum 
to a runaway direction, and focus on the effective potential of the fixed path. 

We choose a path, which is divided into two parts. The first part of the path includes the 
meta-stable vacuum while the other includes a runaway direction. The second path is chosen 
so that it is free of tachyonic modes. It is realized by imposing the condition W^^^Wl = 
on the bosonic mass matrix fl38|l . In other words, it is the condition that the bosonic mass 
matrix coincides with the fermionic one. Then the second part of the path sits on the surface 
determined by the relation W^^^Wl = 0, which is two-dimensional subspace in the seven- 
dimensional space spanned by (0(5), 0(-5), 0(3), 0(-3)7 0(7), 0(-i)! X). The first path is a straight 
line connecting the meta-stable vacuum and a point on the second path. The whole path is 
parametrized by z. The meta-stable vacuum is at ^ = 0, and the runaway direction corresponds 
to 2; — )■ oo. 

The first part is given by 

The second part of the path is given by 

03XoiV3 , /Ms 1 , 03 , / 



Ms ' ^' 27V3A^503^O^ ^ 2iV; 



3<^3 



2 



iV3iV503Xo 3 , fXo 3 X. ^ 2 n / ^ na^ 

' 0(-i) = wiyzr^ ' X = Xoz (1< z), (19) 



M5M7 ' 2M- 



3V3 



where Xq is the VEV of the meta-stable vacuum at the respective temperatures, and Xq = 
Xq{T = 0). 03 is a constantly. 

The effective potential on the path at various temperatures are shown in Fig. [31 The 
vanishing line indicates that tachyonic modes appear on the line (we call this tachyonic line) 
and states can not stay at that line. Therefore the state at the meta-stable vacuum keeps 
staying around the origin {z = 0) and never roll down to the runaway vacuum§|. 

Since X is no longer (pseudo) modulus at states with nonzero VEVs for the mass 
matrices ([6]) and ([7]) are not available to derive the effective potential. The comprehensive 



^ We choose a simply connected path rather than a smoothly connected path like in [Tl] , because even small 
perturbation around ()20p makes tachyonic modes in the practical calculation. 

^ The absolute value of the effective potential at the origin is different from that in Fig. [2j since the size 
of mass matrices increased by two. Due to finite temperature effect, even a massless field contributes to the 
effective potential. This contribution is the constant given by —-^{j^ + ||q) x 2 « —OAIT^. 
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Figure 3: The thermal effective potential on the path. The parameters are M^/M^ = M-j/M^ = 
0.4, N3 = = 1, f = 0.001, and 03 = I/M5. The black circles indicate the value of the origin 
{z = 0) aX X = Xq, = 0. In the interval < 2 < 1 the line is absent if tachyonic modes 
appear. In the plots the numerical value is evaluated with M5 = 1. 

mass matrices Mb and Mp are given in appendix [Bl The above condition W'^^Wl = leads to 

iV30(3)X fM, f 



M5 2iV3iV50{3)^ ' ' 2iV30(3) 



iV3iV50(3)^' ^ fX 



(20) 



M5M7 ' 2M30(3) 

All VEVs in the mass matrices can be taken to be real by field redefinitions without loss of 
generality. One can show that (l20l) are equivalent to (ITTi) and (|T2l) by replacing 0(3) — )• e03 and 
0(5) e~^05- In fact the tree level potential on the surface (120|) is written as 

Taking a limit of 0(3) — 0, X — > 00 and X0(3) — 00 at the same time, Vtree goes to a 
supersymmetric global vacuum. 

The qualitative behavior of the effective potential in Fig. E] is confirmed by an analytic 
calculation without any specification of parameters. At high temperature, the expansion 
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can be applied. The effective potential on the surface fll9p becomes 




01 Mi 



+^ [fMNi + Ni)z' + 0{z-'))+... , 



(22) 



where the Coleman- Weinberg potential is absent as mentioned before. One can see that Vcs 
increases as ~ z"^ for z > 1, and thus the potential is always minimized at 2; = 1 on the second 
part of the path. This suggests that, at high temperature, the stable state exists at least in the 
region z < 1. 

We summarize the evolution of the vacuum structure on the path for the choice of parame- 
ters, iVs = A^5 = 1, M3/M5 = M7/M5 = 0.4, //M| = 0.001 which is considered throughout this 
paper. These parameters allow a vacuum with U{1)r symmetry breaking at zero temperature. 
At high temperature supersymmetry breaking vacuum sits at the origin of the field space {X, z) 
where U{1)r symmetry is not broken. As temperature decreases, potential barrier between the 
origin and the z-direction develops and a local vacuum appears along 2;-direction (Fig. [3]), 
which has a higher potential value than one at the origin. So the origin remains to be favored 
as a vacuum. 

There is a local vacuum along 2;-direction at certain temperatures T = 0.5 ~ 2.0 which is 
separated from the meta-stable vacuum by the tachyonic line as we have mentioned. Mean- 
while, the supersymmetric runaway vacuum develops along z ^ 00 direction around T = 0.5. 
Although the runaway vacuum may be energetically favoured compared to the meta-stable 
vacuum, the local minimum at the origin remains a vacuum state due to the large potential 
barrier between the origin and the runaway vacuum. 

Simultaneously, the local vacuum at the origin of directions (2 = 0) starts to roll down to 
a local vacuum near the origin along the X-direction where the U{1)r symmetry is broken. The 
large potential barrier between X-direction and the directions prevents the f/(l)_R breaking 
vacuum from decaying to the runaway vacuum and keeps it sitting on the meta-stable vacuum. 

5 Conclusions and discussions 

We have studied finite temperature effects on O'Raifeartaigh type models with U{N) global 
symmetry where supersymmetry is spontaneously broken. At zero temperature the model has a 
meta-stable vacuum where U{1)r symmetry is dynamically broken by the one-loop corrections 
for appropriate parameter regions. In addition to a meta-stable vacuum there also exists a 
runaway vacuum. We have investigated the finite temperature effects in the effective potential 
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of the pseudo modulus field X . In order to consider the finite temperature effects along the 
runaway direction, we chose a particular path so that tachyonic modes can be avoided. 
The main features of our study are summarized as follows: 

• Even though we have taken = 1 to study the effective potential, the finite temperature 
effects for the superpotential (ITOi) that allows the global symmetry U{N) is non-trivial. 
The phase structure during the cooling process is highly dependent on the parameter 
matrices M*-^ and N^^ . These matrices can not be chosen freely once one imposes the 
global symmetries. As we have mentioned in the introduction, the global symmetry is a 
key point when one considers the mediation of supersymmetry breaking. 

• Our study reveals the effect of the intermediate and low-temperature region particularly. 
Many of the preceding analyses do not pay much attention to that region, however, when 
one considers the cooling of the universe, it requires rigorous information in the inter- 
mediate and low-temperature region as well as high-temperature. As anticipated from 
the analytic calculations, the U{1)r symmetry is recovered at high temperature which 
is the similar result found in the O'Raifeartaigh model without global symmetries [H]. 
However, we find non-trivial vacuum structures - the landscape - at low temperatures. 
Certain parameter regions that do not have the meta-stable vacua at zero temperature can 
posses U{1)r breaking meta-stable vacuum at finite (low) temperature. This structure is 
not seen in the O'Raifeartaigh model without global symmetries [H]. We obtained this 
result by performing the full numerical analysis without any (high and low temperature) 
approximations . 

• To obtain the correct value of the effective potential in runaway direction, we introduce 
inherently-tachyon-free path where the bosonic mass matrix coincides with the fermionic 
one. Another possible way of the path that is used in [H] is a smooth connection between 
the meta-stable and a runaway vacuum at zero temperature with weighting function. 
That kind of path is not appropriate for our case as pointed out in section 4, since in 
the practical calculation it makes tachyonic modes on the path and causes a meaningless 
analyses. The choice of the path is critical. 

We would like to make some comments about a previous study. In the literature [15], thermal 
effective potential of a class of O'Raifeartaigh model which includes the model of the present 
paper has been analyzed. The author concludes that within certain high-temperature region 
a new meta-stable vacuum that is absent at zero temperature is developed far away from the 
origin X = 0. This kind of vacuum is not considered in section 4. The main reason for this is the 
difference of focusing point and analysis. In [15], the massive scalar fields are integrated out to 
analyze the runaway behaviour X oo, before introducing thermal effects. Then the thermal 
potential is given in analytic form by using high-temperature expansion formula. Meanwhile 
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we keep all the fields explicitly, make the complete mass matrices, calculate eigenvalues of the 
mass matrices, and focus on the effective potential along the fixed path, particularly not far 
from the origin of and X, by using numerical calculation. 

We close this work by pointing out future directions. Here we studied the hidden sector 
that the f/(l)R symmetry breaking is spontaneously broken by the one-loop corrections in 
N" = 1 supersymmetric theory. However, in some JV = 1 supersymmetric model, especially 
accompanying strong gauge dynamics, the perturbative treatment is not available for all regions 
of the pseudo fiat direction. For instance, in the Izawa-Yanagida-Intriligator-Thomas model 
[121 Hn] where an O'Raifeartaigh type sector is dynamically generated by strong gauge dynamics 
in the low energy superpotential, in order to remove degeneracy of pseudo moduli, we have 
to take account of quantum corrections for the Kahler potential. In general, this is a very 
difficult task since the Kahler potential is not holomorphic and thus quantum corrections can 
be estimated at best by perturbative means only in the ultraviolet (weak coupling) region 
of the moduli space parameterizing the pseudo fiat direction which is far from the origin. 
Therefore, the potential behavior in the infrared region remains unclear. However, in an A/^ = 
2 supersymmetric gauge theory one can derive the exact low energy effective action as was 
demonstrated by Seiberg and Witten [201 EI], using the properties of holomorphy and duality. 
In the framework of A/" = 2 supersymmetric gauge theory, some supersymmetry breaking models 
with U{l)ji symmetry breaking were considered [22l |23l [211 (221 126]. It would be interesting to 
study effects of the finite temperature in those models (see also [27]). 
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A Effective potential at high and low temperatures 
A.l High temperature expansion 

In this subsection, we consider series expansions of the thermal functions at high temperature. 
For m/T <C 1, we have the following expansions, 

4 2 2 1 4 2 

JsK/T^] = - + log 



45 12 r2 6 T3 32 ^ a^T^ 

2 \ '+2 

2 I ' (23) 



, /(2/ + l) / 1\ ( 
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and 



J^[m /T ] = - 7T77^ - TTTTT^ log 



360 24 T2 32 T4 ^ afT'^ 



where ( is the Riemann zeta function, = IGvr exp(3/2 — 27^;), a/ = vr exp(3/2 — 27^;) and 
'Je ~ 0.577 is the Euler-Mascheroni constant. Therefore the finite temperature part of the 
effective potential f lT7|) is given by 



vrc, ; 90 V 4 / 24 



■ Ns Np 



127r 

(25) 

where the dots denote terms that are subleading in m/T. A^^, Np are the number of real scalars 
and Weyl fermions and mBi,mpi are mass eigenvalues of the boson and fermion mass matrices. 

In general, it is hard to find the mass eigenvalues analytically. However, one can extract high 
temperature effects without solving eigenvalue equations explicitly. At very high temperature 
T/m ^ 1, the 0(T'^) term is dominant compare to the 0{T) terms in fl25l) □, the coefficient of 

in V^^'^ is given by the trace of mass matrices and m'jp. In our model we find 

+ E^^* = 8(^1 + M| + + NiX' + N^X'), (26) 

i i 

where we have chosen all the parameters real for simplicity. From the expression of the mass 
matrices, it is easy to see that the result is independent of the parameter /. The effective 
potential behaves as 

yW(X, T) ~ T\N^ + Nl)X'^ + 0{T). (27) 

From this expression, one sees that the modulus X tends to be fixed at the origin at high 
temperature if and are not zero. Therefore we expect that the region that allows U{1)r 
breaking vacua is becoming small at high temperature. In fact, this behavior can be observed 
in numerical analysis of the full effective potential (fT7|) (see fig. [1]). 

A. 2 Low temperature expansion 

Let us derive an expression that is useful to study behavior of the effective potential at low 
temperature. For the case m/T ^ 1, the bosonic thermal function is expanded as 



Jb= rdss'f^{-l)^, (2J 

•^0 n=l 



n 



*In (pS)) . 0{T^) term is just the numerical constant and is irrelevant for analyzing the X dependence. 
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where X{s) = e a/'^^+™^^^. The integral over s is rewritten by defining new variable s = 
m/3sinht as 

Jb = E - ( - ^4-) (29) 

I / J \ rlm'^ <n rim I 



n \rfi dm? n dm , 

n=l ^ ' 

This integration is expressed by the modified Bessel functions: 

,3/33 °° 



3 o6 1 

^B = ^ E - iMnm(3) - Ks{nm(3)) , (30) 



4 ^ n 

n=l 

where we have used the expression Kq{z) = J^dt e^^'^°^^^ and relations 

^i^oW = -K,{z), ^,Ko{z) = -\k,{z) - \k,{z). (31) 
Similarly, we have the fermionic thermal function as 

Jf = ^ E iK,{nmP) - K,{nmP)) . (32) 

n=l 

Let us use the asymptotic expansion of the modified Bessel functions 

^1 1 T{iy + l + 1/2) 



2^ ^j2zynT{u -1+1/2) 

^ 1 e-^fl + ^4^ + ---V (kl^oo), (33) 



2 V 8^ 

where the dots are subleading order in z~K Then, we have 

- -y|(m/3)iLi|(e-^). (34) 

Here Lis{z) = YlhLi polylogarism function and we have used the fact that the dominant 

contribution comes from n = 1 sector in the expansions (1301) and (152]) . Again, using the 
asymptotic behavior Lis(2;) = z + 2~^z'^ + ?)~'^z^ + ■ ■ ■ {\z\ -+ 0), we find the low-temperature 
expansion of the thermal functions 

TT , 3 



JB{mP) ^|(^/3)i(^e-'"/5 + 2"5e-2"^^ + 3-^6-3™'^ + (35) 

Jp^mP) ~ -y|(m/3)i (-e^"'^ + 2-ie-2"^'^-3-te-3"^'' + ---) . (36) 

As we have discussed, at the origin of </>*, the mass for cff is given by m ~ M + NX. Therefore 
we expect that the finite temperature part of the effective potential at low temperature behaves 
like 

\/«(X,T) ~ -T^Tr |((M + iVX)T-i)ie-(*''+^^)^"'} . (37) 
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This expression would be useful for a calculable model to see the thermal evolution of vacuum 
analytically at low temperature. In our case, since analytic forms of eigenvalues are complicated, 
it is difficult to see the evolution analytically. However, the low temperature expansion is still 
useful for numerical calculations. It needs shorter time to evaluate the effective potential 
compare to the full thermal functions f llSp . fll6p . One can use the low temperature expansion 
to check the full result in the numerical analysis in section 3. 



B Mass matrices 

The general mass matrices are given as follows. 



Mi 



Mi 



wlw^^ 

W'^Wl- 



(38) 
(39) 



where 



I 



M|+7V|{|Xp 

+ I0-5P) 













hhN-i<f>3+NlX^r, 




A/2+Ar2(|x|2 




NsiJVhX"! 

+N3<l>3<t>l) 








7V20_5Xt 


+N5<j)-5<pl3) 




A/|+iV|(|X|2 

+ I<^-3P) 


Ni<t>3<t>U 
















M|+7V|(|X|2 





M-iNzX'i 


A/5Ar5</._5+A'|Xt</,„3 


M7N5X 











A/2 
















M3N3X 





Ml 


M3N3<j>-3 






NiX4>l 








JV|(|</'3P + |05P 
+ |</._3|2 + |</,_5|2) 



(40) 



15 






N5{f+N3<l,3<P-3 














NsiNsXcPs 














+M7<t>7) 


Ns{f+N3<l^3'P-3 

















N5{N5X<I>.5 














+M50_3) 



















N3iN3X<j,3 














+M5<t>5) 








N-i{f+N3<l>3'l>-3 











N3iN3X<f>_3 






+N5<t>5<t>-5) 








+M3<f>-l) 
















































N3{N3X(I>3 













+M^<|>^) 


+h'h<l,-3) 




+M3</>-l) 









(41) 

X and (j)r are VEVs of X and 0(r.) respectively. 
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